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We consider the possibility of several different mechanisms contributing to the (/3/3)o.y-decay 
amplitude in the general case of CP nonconservation: light Majorana neutrino exchange, heavy 
left-handed (LH) and heavy right-handed (RH) Majorana neutrino exchanges, lepton charge non- 
conserving couplings in SUSY theories with i?-parity breaking. If the (/?/3)oi/-decay is induced by, 
e.g., two "non- interfering" mechanisms (light Majorana neutrino and heavy RH Majorana neutrino 
exchanges), one can determine \r]if and \r]jf, rji and rjj being the two fundamental parameters 
, characterising these mechanisms, from data on the half-lives of two nuclear isotopes. In the case 

when two "interfering" mechanisms are responsible for the (/3/3)o,y-decay, [rjif and \r}jf and the 



interference term can be uniquely determined, in principle, from data on the half-lives of three nu- 
clei. Given the half-life of one isotope, the "positivity conditions" \rii\^ > and > lead to 

(-H . stringent constraints on the half- lives of the other (/?/3)oi/-decaying isotopes. These conditions, as 

well as the conditions for constructive (destructive) interference are derived and their implications 
are analysed in two specific cases. The experimental limits on neutrino masses obtained in the '^H 
/3-decay experiments can constrain further the multiple mechanisms of (/3/3)oi/-decay if one of the 
^ ' mechanisms involved is the light Majorana neutrino exchange. The measurements of the half-lives 

\l ' with rather high precision and the knowledge of the relevant nuclear matrix elements with rela- 

tively small uncertainties is crucial for establishing that more than one mechanisms are operative 
in (/?/3)oiy-decay. The method considered by us can be generalised to the case of more than two 
(/9/3)oi/-decay mechanisms. It allows to treat the cases of CP conserving and CP nonconserving 
fT^ , couplings generating the (/?/3)oi/-decay in a unique way. 

o 
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I. INTRODUCTION 
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X 



If neutrinoless double beta ((/J/3)oi/-) decay will be observed, it will be of fundamental importance 
to determine the mechanism which induces the decay. We know that neutrinos have mass and mix, and 
■ if they are Majorana particles they should trigger the decay at some probability level. The fundamental 

parameter which controls the (/3/3)oi/-decay rate in this case is the effective Majorana mass: 

light 

<m> = {Uejfmj , (all mj > 0) , (1.1) 

j 

where U is the Pontecorvo, Maki, Nakagawa and Sakata (PMNS) neutrino mixing matrix [ll-[l| and rrij 
are the light Majorana neutrino masses, mj ^ 1 eV. The (/3/3)oi/-decay rate depends on the type of 
neutrino mass spectrum which can be hierarchical, with partial hierarchy or quasi-degenerate (see, e.g., 
Q). Using the data on the neutrino oscillation parameters it is possible to show (see, e.g., |5i]) that in the 
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case of normal hierarchical spectrum one has |<m>| ^ 0.005 eV, while if the spectrum is with inverted 
hierarchy, 0.01 eV !^ \<m>\ ^ 0.05 eV. A larger value of |<m>| is possible if the light neutrino mass 
spectrum is with partial hierarchy or of quasi-degenerate type. In the latter case |<to>| can be close 
to the existing upper limits. 

The most stringent upper hmits on |<to>| were set by the IGEX Q, CUORICINO and NEMOS 
experiments with ^^Ge, ^^°Te and ^"^^Mo, respectively ^. The IGEX collaboration has obtained 
for the half-life of ^^Qe, T°J^ > 1.57 x lO^^ yr (90% C.L.), from which the limit |<m>| < (0.33 - 
1.35) eV was derived Using the recent more advanced calculations of the corresponding nuclear 
matrix elements (including the relevant uncertainties) fl^l one finds: |<m>| < (0.22 — 0.35) eV. The 
NEM03 and CUORICINO experiments, designed to reach a sensitivity to |<m>| (0.2 - 0.3) eV, set 
the limits: |<m>| < (0.61- 1.26) eV [8] and |<m>| < (0.19 - 0.68) eV (90% C.L.), where estimated 
uncertainties in the NME are accounted for. The two upper limits were derived from the experimental 
lower limits on the half-lives of i°°Mo and ^^Oxe, T^J^ > 5.8 x 10^3 yr (90%C.L.) [8] and T^J^ > 

3.0 X 10^4 yr (90%C.L.) 0. With the NMEs and their uncertainties calculated in jl^, the NEM03 and 
CUORICINO upper limits read, respectively: |<m>| < (0.50-0.96) eV and |<to>| < (0.25-0.43) eV. 
The best lower Hmit on the half-Hfe of ^'^Ge, T^J^ > 1.9 x 10^^ yr (90% C.L.), was found in the Heidelberg- 
Moscow ^^Ge experiment (H. It corresponds to the upper limit [l3 |<m>| < (0.20 — 0.35) eV. A positive 
(/3/3)o,.-decay signal at > 3a, corresponding to T^J^ = (0.69-4.18) x 10^^ yr (99.73% C.L.) and implying 
|<m>| = (0.1 — 0.9) eV, is claimed to have been observed in [lo|, while a later analysis reports evidence 
for (/3/3)oi/-decay at 6cr corresponding to |<m>| = 0.32 ± 0.03 eV [ll|. 

Most importantly, a large number of projects aim at a sensitivity to |<to>| ^ (0.01 — 0.05) eV fll|: 
CUORE (i30Te), GERDA (^^Ge), SuperNEMO, EXO C^^Xe), MAJORANA (^^Ge), MOON (lo^Mo), 
COBRA ("^Cd), XMASS ("*^Xe), CANDLES (^**Ca), KamLAND-Zen (^^exe), SNO+ {^^°Nd), etc. 
These experiments, in particular, will test the positive result claimed in fTTj. Let us note that the 
measurement of |<m>| can provide unique information on the absolute scale of neutrino masses, the 
type of neutrino mass spectrum and the Majorana phases in the PMNS matrix [14J . 

The light Majorana neutrino exchange can be called the "standard" mechanism of the {/3(3)o^-deca,y. 
The observation of (/3/3)oj/-decay would imply that the total lepton charge L is not conserved. This 
would also imply that the massive neutrinos get a Majorana mass [15] and therefore are Majorana 
particles (see, e.g. [l^). However, the latter does not guarantee that the dominant mechanism inducing 
the (/3/3)oi/-decay is the light Majorana neutrino exchange since the Majorana mass thus generated is 
exceedingly small. The (/3/3)o;y-decay can well be due to the existence of interactions which do not 
conserve the total lepton charge L, AL = ±2. A number of such interactions have been proposed in 
the literature: heavy Majorana neutrinos coupled to the electron in the V — A charged current weak 
interaction Lagrangian, supersymmetric (SUSY) theories with i?-parity breaking terms which do not 
conserve the total lepton charge L, L-nonconserving couplings in the Left-Right symmetric theories, etc. 
At present we do not have evidence for the existence of AL ^ terms in the Lagrangian describing the 
particle interactions. Nevertheless, such terms can exist and they can be operative in the (/3/3)oi/-decay. 
Moreover, it is impossible to exclude the hypothesis that, if observed, the (/3/3)oi/-decay is triggered by 
more than one competing mechanisms. 

The possibility of several different mechanisms contributing to the (/3/3)ojy-decay amplitude was 
considered recently in (l3| assuming that the corresponding AL — ±2 couplings are CP conserving. 
By exploiting the dependence of the nuclear matrix elements on the decaying nucleus, it was shown 
that, given the experimental observation of the (/3/3)oi/-decay of sufficient number of nuclei, one can 
determine and/or sufficiently constrain the fundamental parameters associated with the lepton charge 
nonconserving couplings generating the (/3/3)oiy-decay. 

The present work is a natural continuation of the study performed in [I^ . We consider the possibility 
of several different mechanisms contributing to the (/3/3)o[/-decay amplitude in the general case of CP 
nonconservation: light Majorana neutrino exchange, heavy left-handed (LH) and heavy right-handed 
(RH) Majorana neutrino exchanges, lepton charge non-conserving couplings in SUSY theories with im- 
parity breaking. If the (/3/3)oi/-decay is induced by, e.g., two "non-interfering" mechanisms (light Majorana 
neutrino and heavy RH Majorana neutrino exchanges), one can determine the absolute values of the two 
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fundamental parameters, characterising these mechanisms, from data on the half-Hves of two nuclear 
isotopes. In the case when two "interfering" mechanisms are responsible for the (/3/3)o;y-decay, the absolute 
values of the two relevant parameters and the interference term can be uniquely determined from data 
on the half-lives of three nuclei. In the specific examples considered of two "noninterfering" and two 
"interfering" mechanisms, namely, the light Majorana neutrino and the heavy RH Majorana neutrino 
exchanges, and the light Majorana neutrino and the dominant gluino exchanges, we present illustrative 
examples of determination of the relevant fundamental parameters and of possible tests of the hypothesis 
that more than one mechanism is responsible for the (/3/3)o;y -decay, using as input hypothetical half-lives 
of '^^Ge, ^^'^Te and ^""^Mo. The effects of the uncertainties in the values of the nuclear matrix elements 
(NMEs) on the results of the indicated analyzes are also discussed and illustrated. 

The method considered by us can be generalised to the case of more than two (/3/3)oi/-decay mecha- 
nisms. It has also the advantage that it allows to treat the cases of CP conserving and CP nonconserving 
couplings generating the {/3/3)oi/-decay in a unique way. 



II. DIFFERENT MECHANISMS OF (/3^)o^-DECAY 

We will consider in the present article the following mechanisms of (/3/3)o;y-decay: the exchange of 
light Majorana neutrinos; the exchange of heavy "left-handed" (LH) Majorana neutrinos; the exchange of 
heavy "right-handed" (RH) Majorana neutrinos; and two mechanisms associated with possible i?-parity 
breaking in SUSY theories. Below we discuss briefly the lepton number violating (LNV) parameters and 
the nuclear matrix elements associated with each of the indicated mechanisms. 

Assuming the dominance of a single LNV mechanism characterised by a parameter 7?^ , where the 
index k denotes the mechanism, the inverse value of the (/3/?)o;y-decay half-life for a given isotope {A, Z) 
can be written as 

^ = \^^^^y\^G''^{E,^Z)\M":\\ (2.1) 

^1/2 

where G'^'^{Eo,Z) and M'^" are, respectively, the known phase-space factor {Eq is the energy release) 
and the nuclear matrix element of the decay. The latter depends on the mechanism generating the decay 
and on the nuclear structure of the specific isotopes {A, Z), {A, Z + 1) and {A, Z + 2) under study. 

The phase space factors G'^'^{Eo, Z), which include the fourth power of the "standard" value of the 
axial-coupling constant — 1.25, are tabulated in ref. |18|]; for ^^Ge, ^'^Se, ^'^"Afo and ^'^"Te are given 
in Table I. For a given isotope (A, Z), G^'^{Eo, Z) contains the inverse square of the nuclear radius R{A) 
of the isotope, R~^{A), compensated by the factor R{A) in M' ^ . The assumed value of the nuclear 
radius is R[A) = tqA^^^ with tq = 1.1 fm. 

The nuclear matrix element M' ^ is defined as 

M":=(^)\llr (2.2) 

This definition of Af'^'' [IqI allows to display the effects of uncertainties in qa and to use the same phase 
factor G'^''{Eo, Z) when calculating the (/?/3)o,y-decay rate. 



A. Light Majorana Neutrino Exchange 

In the case of the light Majorana neutrino exchange mechanism of (/?/3)oi/-decay, the LNV parameter 
is given by: 

<m> , ^ 

Vu = , 2.3 

where <m> is the effective Majorana mass (see, e.g., flG^). Under the assumption of n light massive 
Majorana neutrinos coupled to the electron in the weak charged lepton current, the effective Majorana 
mass is given in eq. (|l.ip . Thus, <to> depends on the elements of first row of the PMNS neutrino 



mixing matrix, Uej, j — 1, 2.3, ... The PMNS matrix U is not assumed to be CP conserving and at least 
two of the elements Uej contain physical CP violating phases pol W\\ (see also, e.g., @). In the case of 
3 light neutrinos and the "standard" parametrisation of C/ jsl , the elements Ue2 and Ues contain the two 
physical CP violating Majorana phases [20] and Ue3 contains the Dirac phase as well. 

The expression for <m>, eq. corresponds to the contribution from the standard {V — A) 

charged current (CC) weak interaction. The nuclear matrix element M^'^ for different isotopes (A, Z) is 
given in [Ullll (see also Table I). 



B. Heavy Majorana Neutrino Exchange Mechanisms 



We assume that the neutrino mass spectrum includes, in addition to the three light Majorana 
neutrinos, heavy Majorana states Nk with masses Mk much larger than the typical energy scale of the 
(/3/?)oi/-decay, Mk S> 100 MeV; we will consider the case of Mk ^ 10 GeV. Such a possibility arises if the 
weak interaction Lagrangian includes right-handed (RH) sterile neutrino fields which couple to the LH 
flavour neutrino flelds via the neutrino Yukawa coupling and possess a Majorana mass term. The heavy 
Majorana neutrinos Nk can mediate the (/3/3)o,y-decay j23j similar to the light Majorana neutrinos via the 
V — A charged current weak interaction. The difference between the two mechanisms is that, unlike the 
light Majorana neutrino exchange which leads to a long range inter-nucleon interactions, in the case of 
Mk ^ 10 GeV of interest the momentum dependence of the heavy Majorana neutrino propagators can be 
neglected (i.e., the Nk propagators can be contracted to points) and, as a consequence, the corresponding 
effective nucleon transition operators are local. The LNV parameter in the case when the (/3/?)oi/-decay 
is generated by the {V — A) CC weak interaction due to the exchange of Nk can be written as: 
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(2.4) 



where nip is the proton mass and Uek is the element of the neutrino mixing matrix through which Nk 



V -A 



V + A 



W, 



N 



kR 



Wr 



V -A 



V + A 



Figure 1: Feynman diagrams for the (/3/3)oi^-decay, generated by the light and heavy LH Majorana neutrino 
exchange (left panel) and the heavy (RH) Majorana neutrino exchange (right panel). 

couples to the electron in the weak charged lepton current. We note that jry^ | is suppressed by both the 
ratio TOp/Mfc and the magnitude of Uek (see, e.g., [23|). 

If the weak interaction Lagrangian contains also {V + A) (i.e., right-handed (RH)) charged currents 
coupled to a RH charged weak boson Wr, as, e.g., in the Left-Right Symmetric theories, we can have 
also a contribution to the (/?/3)oi/-decay amplitude generated by the exchange of virtual Nk coupled to 
the electron in the hypothetical (V + A) CC part of the weak interaction Lagrangian [2^. In this case 
the corresponding LNV parameter can be written as: 




Here Vek are the elements of a mixing matrix by which Nk couple to the electron in the {V + A) charged 
lepton current ^, Mw is the mass of the Standard Model charged weak boson, Mw — 80 GeV, and Mwr 
is the mass of Wr. It follows from the existing data that Mwr ^ 2.5 TeV. Thus, |?7^| is suppressed 
by the factor {Mw/MwRf- 

If CP invariance does not hold, which we will assume to be the case in what follows, Uek and Vek 
will contain physical CP violating phases at least for some k and thus the parameters j]^ and rj^ will 
not be real. 

As can be shown, the nuclear matrix elements corresponding to the two mechanisms of (/3/3)o;y-decay 
with exchange of heavy Majorana neutrinos Nk, described in the present subsection, are the same and 
are given m [18]. We will denote them by M^"" (and M'"^). 

Finally, it is important to note that the current factor in the (/3/3)oi/-decay amplitude describing the 
two final state electrons, has different forms in the cases of (/3/3)o,y-decay mediated by {V — A) and by 
{V + A) CC weak interactions, namely, e(l + 75)e'^ = 2eL ie'^)R and e(l — 75)6'^ = 2e~R (e'^)^, respectively, 
where e'^ = C(e)^, C being the charge conjugation matrix (see, e.g., [16]). The difference in the chiral 
structure of the two currents leads to a specific phase space factor of the interference term in the rate of 
(/3/3)oj/-decay, triggered by two mechanisms whose respective contributions to the (/3/3)o;/-decay amplitude 
involve the two different electron current factors. The phase space factor of the interference term under 
discussion is significantly smaller than the phase space factors of the contributions to the (/3/3)oi/-decay 
rate due to each of the two mechanisms, which leads to a relatively strong suppression of the interference 
term [23] (see further). 

C. SUSY Models with R-Parity Non-conservation 

The SUSY models with R-parity non-conservation include LNV couplings which can trigger the 
(/3/?)oi/-decay. Let us recall that the R-parity is a multiplicative quantum number defined by i? = 
(^_l^2S+3B+L ^ where S, B and L are the spin, the baryon and lepton numbers of a given particle. The 
ordinary (Standard Model) particles have i? = +1, while their superpartners carry 




Figure 2: Feynman diagrams for (/3/3)oi^-decay due to the gluino exchange mechanism. 

i? = — 1. The LNV couplings emerge in this class of SUSY models from the R-parity breaking (^p) part 
of the superpotential 

Wfi^ = \,kUL,El + K^^UQ.Dl + ^i,UH2, (2.6) 



^ We have neglected the contributions to and, more generally, to the (/3/3)oi/-decay amplitude due to the possible but 
small mixing between W and Wr bosons. 



where L, Q stand for lepton and quark SU{2)l doublet left-handed superfields, while E'^, D'^ for lepton 
and down quark singlet superfields. Here we concentrate only on the trilinear A'-couplings. The A'- 
couplings of the first family of particles and sparticles relevant for (/3/3)o,y-decay are given in terms of the 
fields of the LH electron, electron neutrino UeL, LH selectron and sneutrino VeL, LH and RH u- and 
c?-quarks, u^ ji and d^.i?, and LH and RH u- and c?-squarks, ml,_r, (Ilm, by: 
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(2.7) 



At the quark-level there are basically two types of SUSY mechanisms of (/?/3)o^-decay: a short- 
range one with exchange of heavy Majorana and scalar SUSY particles (gluinos and squarks, and/or 
neutralinos and selectrons) [28-33], and a long-range mechanism involving the exchange of both heavy 
squarks and light Majorana neutrinos (sil - lssj . We will call the latter the "squark-neutrino" mechanism. 



The Case of Gluino Exchange Dominance 

Assuming the dominance of the gluino exchange in the short-range mechanism, one obtains the 
following expression for the corresponding LNV parameter: 



Af 
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(2.8) 



Here, Gp is the Fermi constant, — g|/(47r), 173 being the SU(3)c gauge coupling constant, m^^, ra^^ 
and nig are masses of the LH u-squark, RH d-squark and gluino, respectively. 

The nuclear matrix element associated with the gluino exchange mechanism, M°r, was calculated 
in (39ll40|. The electron current factor in the term of the (/3/3)o!y-decay amplitude corresponding to the 
gluino exchange mechanism under discussion has the form e(l -1-75)6'^ = 2eL (e'^)jj, i.e., it coincides with 
that of the light (or heavy LH) Majorana neutrino exchange. Thus, when calculating the (/3/3)ojy-decay 
rate, the interference between the two terms present in the (/3/3)ojy-decay amplitude, corresponding to 
the indicated two mechanisms, has the same phase space factor as the contributions due to each of the 
two mechanisms. As a consequence, the interference term will not suffer from phase space suppression. 



The S quark- Neutrino Mechanism 

In the case of squark-neutrino exchange (ssj . the (/3/3)oi/-decay amplitude does not vanish in the limit 
of zero Majorana neutrino mass, in contrast to the case of the "standard" light (LH) Majorana neutrino 
exchange. This is a consequence of the chiral structure of the corresponding I^p SUSY couplings. The 
contribution due to the squark-neutrino exchange to the (/3/3)oi/-decay amplitude is roughly proportional 
to the momentum of the virtual neutrino, which is of the order of the Fermi momentum of the nucleons 
inside of nucleus, pp ~ 100 MeV. The corresponding LNV parameter is defined as |38| 



- Ev^-n20f,. -i . (2.9) 



Here we use the notations d(^k) = d, s,b and assumed that there are 3 light Majorana neutrinos. This 
LNV parameter vanishes in the absence of dkL — dkR - mixing, i.e., when 6'^ = 0. The nuclear matrix 
element for the squark-neutrino mechanism, M?", is given in [38|. 



III. NUCLEAR STRUCTURE CALCULATIONS 

In what follows the (/3/3)oi/-decay nuclear matrix elements M'^J , -Wjv i M''^y and M'^^ are evaluated 
for ^^Ge, ^"^Se, ^^^Mo and ^^^Te. These nuclei are considered as most probable candidate sources for 
the next generation of the experiments searching for (/3^)oi/-decay. 



Table I: The phase-space factor G^" {Eq, Z) and the nuclear matrix elements M'^J (light Majorana neutrino 
exchange mechanism), M' ^ (heavy Majorana neutrino exchange mechanism), M' y (mechanism of gluino ex- 
change dominance in SUSY with trilinear R-parity breaking term) and M' ^ (squark- neutrino mechanism) for 
the (/3/3)o;.-decays of '^^Ge, Se, ^°°Mo and ^^°Te. The nuclear matrix elements were obtained within the 
Self-consistent Renormalized Quasiparticle Random Phase Approximation (SRQRPA). See text for details. 
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We used the Self-consistent Renormalized Quasiparticle Random Phase Approximation 
(SRQRPA) |4l| to calculate the nuclear matrix elements (NMEs) of interest. The SRQRPA takes into 
account the Pauli exclusion principle and conserves the mean particle number in correlated ground state. 

For each of the four nuclei, two choices of single-particle basis are considered. The intermediate size 
model space has 12 levels (oscillator shells N=2-4) for ^^Ge and *^S'e, 16 levels (oscillator shells N=2-4 
plus the f+h orbits from N=5) for ^°'^A/o and 18 levels (oscillator shells N=3,4 plus f+h+p orbits from 
N=5) for ^^^Te. The large size single particle space contains 21 levels (oscillator shells N=0-5) for ^^Ge, 
^^Se and -'^™Mo, and 23 levels for ^^'■^Te (N=l-5 and i orbits from N=6). In comparison with previous 
studies we omitted the small space model which is not sufficient to describe realistically the tensor 
part of the (/3/?)oiy-decay nuclear matrix elements. 

The single particle energies were obtained by using a Coulomb-corrected Woods-Saxon potential. 
Two-body G-matrix elements we derived from the Argonne and the Charge Dependent Bonn (CD-Bonn) 
one-boson exchange potential within the Brueckner theory. The schematic pairing interactions have 
been adjusted to fit the empirical pairing gaps j42|. The particle-particle and particle-hole channels of 
the G-matrix interaction of the nuclear Hamiltonian H are renormalized by introducing the parameters 
Qpp and Qph, respectively. The calculations have been carried out for gph = 1.0. The particle-particle 
strength parameter gpp of the SRQRPA is fixed by the data on the two-neutrino double beta decays 
(lOl [22|. In the calculation of the (/?/3)oi/-decay NMEs. the two-nucleon short-range correlations derived 
from same potential as residual interactions, namely from the Argonne or CD-Bonn potentials, were 
considered |43j . 

The calculated NMEs M'l" , M'°^ , M'°y and M'f are listed in Table HI We see that a significant 
source of uncertainty is the value of the axial- vector coupling constant qa- Further, the NMEs associated 
with heavy neutrino exchange are sensitive also to the choice of the NN interaction, the CD-Bonn or 
Argonne potential. These types of realistic NN interaction differ mostly by the description of the short- 
range interactions. 

Finally, we notice that all NMEs given in Table |T] are real and positive. 



IV. ANALYSIS 



We illustrate the possibility to get information about the different LNV parameters when two or 
more mechanisms are operative in (/3/3)oi/-decay, analysing the following two cases. First we consider two 
competitive "not-interfering" mechanisms of (/3/3)o^-decay: light left-handed Majorana neutrino exchange 
and heavy right-handed Majorana neutrino exchange. In this case the interference term arising in the 
(/3/3)o;y-decay half-life from the product of the contributions due to the two mechanisms in the (/3/3)oi/- 
decay amplitude, is strongly suppressed [27| as a consequence of the different chiral structure of the final 
state electron current in the two amplitudes. The latter leads to a different phase-space factor for the 
interference term, which is typically by a factor of 10 smaller than the standard one (corresponding to 
the contribution to the (/3/3)o;y-decay half-life of each of the two mechanisms). More specifically, the 
suppression factors for "^^Ge, ^'^Se, i°°Mo and ^^°Te read, respectively (l^: 0.13; 0.08; 0.075 and 0.10. 
It is particularly small for ''^ Co: 0.04. In the analysis which follows we will neglect the contribution of 
the interference term in the (/3/3)o;y-decay half-life. The effect of taking into account the interference term 
on the results thus obtained, as our numerical calculations have shown, does not exceed approximately 
10%. 

In the case of negligible interference term, the inverse value of the (/3/3)oiy-decay half-life for a given 
isotope (A,Z) is given by: 



r^A'WZ? + \mm%\\ (4.1) 



where the index i denotes the isotope. The values of the phase space factor G^{E, Z) and of the NMEs 
M'^'^^ and M'°^j^ for ^^Qe, ^'^Se, ^°°Mo and "°Te are Hsted in TableHl The parameters \-q^\ and \'qR\ are 
defined in eqs. and (P3|l . 

In the second illustrative case we consider (/3/3)o,y-decay triggered by two active and "interfering" 
mechanisms: the light Majorana neutrino exchange and the gluino exchange. In this case, for a given 
nucleus, the inverse of the (/3/3)o,y-decay half-life is given by: 



T0;2_,G^(i?,Z) 



M''\M'Z? + I^Af iMtvP + 2cosa\M%\\M'tM\m'\- (4.2) 



Here \r}xi \ is the basic parameter of the gluino exchange mechanism defined in eq. (j2.8p and a is the 
relative phase of rjxi and rjy . The values of the NMEs of the mechanisms considered are listed in Table U 
In the illustrative examples of how one can extract information about \'qu\i etc. we use as input 
hypothetical values of the (/3/3)oy-decay half-life of ^^Ge satisfying the existing lower limits and the value 
claimed in ref. fll 01, as well as the following hypothetical ranges for T^f^{^°°Mo) and T'j');2(^^°Te): 



T°;2(^^Ge) > 1.9 X lO^Sy, r{';2(^^Ge) = 2.23^°;;^^ x 10'^2/ 
5.8 X lO^^y < T^J^C^°Mo) < 5.8 x 10^^, 3.0 x lO^'^y < T^J^C^^Te) < 3.0 x lO^^y 



(4.3) 



Let us note that 5.8 x lO'^'^ y and 3.0 x 10^'* y are the existing lower bounds on the half-lives of ^"^Afo 
and i30Te Hi. 

In the analysis which follows we will present numerical results first for gA — 1-25 and using the 
NMEs calculated with the large size single particle basis ("large basis") and the Charge Dependent Bonn 
(CD-Bonn) potential. Later results for gA — 1.0, as well as for NMEs calculated with the Argonne 
potential, will also be reported. 

As we will see, in certain cases of at least one more mechanism being operative in (/3/3)o,y-decay 
beyond the light neutrino exchange, one has to take into account the upper limit on the absolute scale of 
neutrino masses set by the '^H /?-decay experiments [44. 145;] : m{De) < 2.3 eV. In the case of (/3/3)o;y-decay, 
this limit implies a similar limit on the effective Majorana mass | < m > | < 2.3 eV. The latter translates 



We remind the reader that for mi, 2, 3 ^ 0.1 eV the neutrino mass spectrum is quasi-degenerate (QD), mi = m,2 = mg i 



>> Am|j^, |Amgj|. In this case we have m(uc) = m and |<m>| ^ m. 



into the following limit on the conveniently rescaled parameter |77i/p: 

X 10^° < 0.21. (4.4) 

A more stringent limit on the absolute neutrino mass scale and therefore on | < m > | is planned to 
be obtained in the KATRIN experiment jUj: |<to>| < 0.2 eV (90% C.L.). This corresponds to the 
following prospective limit on \ri^\'^: 

Iry^P X 10^° < 1.6 X 10"^ (4.5) 

As the results presented in Section 6 indicate, if the limit of 0.2 eV will be reached in the KATRIN 
experiment, this will lead to severe constraints on some of the solutions for |?7y ^ obtained in the case of 
two "interfering" mechanisms, one of which is the light neutrino exchange. 



V. TWO "NON-INTERFERING" MECHANISMS 



In the case of two active non-interfering (/3/3)o;y-decay generating mechanisms, which we will assume 
to be the light LH and heavy RH Majorana neutrino exchanges [131 , it is possible to extract the absolute 
values of the corresponding two LNV fundamental parameters Irji^l and \iir\, using the "data" on the half- 
lives of two different nuclei undergoing the (/3^)oj/-decay. Indeed, using eq. (|4.ip we can set a system of 
two linear equations for two unknowns using as input the two half-lives: 

= KmZ\' + \vBm%\', = m'\mZ? + immtN?- (5.1) 

The solutions read: 

^ \M%\yT,G, - \M%\yT,G2 \M%\^/T^G^-\M%\'/T,G, 



2,1^1 



Obviously, solutions giving \rj^\'^ < and/or \r]fj\'^ < are unphysical. Given a pair of nuclei {Ai,Zi), 
{A2, Z2) of the three ^^Ge, ^"^Mo and ^'^"Te we will be considering, and Ti, and choosing (for convenience) 
always Ai < A2, positive solutions for and |?7_rP are possible for the following range of values of T2: 

T,G,\M%\^ ^ TiGilMr.P 

n < I2 < TT- , 

G2\M%\^ - G2|M'°;;|2 

where we have used the fact that |M'°|;| V|M'2^^|2 > {M'l'^j^l^ /{M'^'^j^l^ (see TableU) ^. Using the values 
of the phase-space factors and the two relevant NME given in Table I in the columns "CD-Bonn, large, 
gA = 1.25", we get the positivity conditions for the ratio of the half- lives of the different pairs of the 
three nuclei of interest: 

T0-,(1"0Mo) TO-,("Ore) T^YJ^'^'Te) 

^ %(W - ' - - ' - - ■ ^'-'^ 

In the case of gA ~ 1.0 we find: 

rOr,(i°°Mo) T0-(i30Te) T^^J^^OTe) 

0-15 <^4rp^< 017, 0.17 <^^^^^ (5.5) 



This condition will exhibit a relatively weak dependence on the value of gA in the cases of mechanisms in which the 
Gamow- Teller term in the NMEs dominates (as in the gluino exchange dominance and the squark-neutrino exchange 
mechanisms). Indeed, the factor (1.25)* (corresponding to = 1-25) and included in the definition of the phase space 
terms Gi,2, cancels in the ratio G1/G2, and | Af'i_^(iv) | V|M'2,^(]V) P = I A^i°t(JV) I Vl ^^^.'^(jv) 1^ (see eq. 



It is quite remarkable that the physical solutions are possible only if the ratio of the half-lives of all the 
pairs of the three isotopes considered take values in very narrow intervals. This result is a consequence 
of the values of the phase space factors and of the NME for the two mechanisms considered. In the case 
of the Argonne potential, "large basis" and gA = 1-25 (1.0) (see Table |T|) we get very similar results: 

0.15 < ^4rp^ < 0.18 , 0.18 < .^^^ < 0.24 (0.25) , 1.22 < ^^^p^ < 06 (1.42) . 

(5.6) 

If it is experimentally established that any of the three ratios of half-lives considered lies outside the 
interval of physical solutions of \ri^\'^ and |??_rP, obtained taking into account all relevant uncertainties, one 
would be led to conclude that the (/3/3)oi/-decay is not generated by the two mechanisms under discussion. 
In order to show that the constraints given above are indeed satisfied, the relevant ratios of (/3/3)o,y-decay 
half-lives should be known with a remarkably small uncertainty (not exceeding approximately 5% of the 
central values of the intervals). 

Obviously, given the half- life of one isotope, constraints similar to those described above can be 
derived on the half- life of any other isotope beyond those considered by us. Similar constraints can be 
obtained in all cases of two "non-interfering" mechanisms generating the (/3/3)o;/-decay. The predicted in- 
tervals of half- lives of the various isotopes will differ, in general, for the different pairs of "non-interfering" 
mechanisms. However, these differences in the cases of the of the (/3/3)o;/-decay triggered by the exchange 
of heavy Majorana neutrinos coupled to (V+A) currents and i) the gluino exchange mechanism, or ii) 
the squark-neutrino exchange mechanism, are extremely small. One of the consequences of this feature 
of the different pairs of "non-interfering" mechanisms considered by us is that if it will be possible to rule 
out one of them as the cause of (/3/3)o;y-decay, most likely one will be able to rule out all three of them. 
The set of constraints under discussion will not be valid, in general, if the (/3/3)oi/-decay is triggered by 
two "interfering" mechanisms with a non-negligible interference term, or by more than two mechanisms 
with significant contributions to the (/3/3)oi/-decay rates of the different nuclei. 

Evidently, if one of the two solutions is zero, for example |77_r,P = 0, then only one of the two (/3/3)oiy- 
decay mechanisms is active. Since for the two mechanisms considered we have (M' i'^^^)'^ (M' 2'^^)'^ — 
{M'^YiM'^'^^)^ 7^ 0, the condition that \riR\^ = reads: 

\M%\^ Ti Gi = \M'Z\' T2 G2 , \vr\^ = . (5.7) 
The condition that \r]^\'^ = has a similar form: 

\M'%\' Ti Gi = \M'fj,\^ T2 G2 , = . (5.8) 



If only the light neutrino exchange mechanism is present and the NME are correctly calculated, the 
(/3/3)oi/-decay effective Majorana mass (and |?7;Jf ) extracted from all three (or any number of) (/3/3)oi/- 
decay isotopes must be the same (see, e.g., [46], S^). Similarly, if the heavy RH Majorana neutrino 
exchange gives the dominant contribution, the extracted value |?7_rP must be the same for all three (or 
more) (/3/3)oi,-decay nuclei. 

We analyze next the possible solutions for different combinations of the half-lives of the following 
isotopes: ^^Ge, ^"^Mo and ^■^°Te. Assuming the half-lives of two isotopes to be known and using the 
physical solutions for \ri^\'^ and |?7_rP obtained using these half-lives, one can obtain a prediction for the 
half- life of the third isotope. The predicted half-life should satisfy the existing lower limits on it. In the 
calculations the results of which are reported here, we fixed the half-life of one of the two isotopes and 
assumed the second half- life lies in an interval compatible with the existing constraints. We used the 
value of T0^2(^^G'e) ^nd values of ^^''/^(^^^ Mo) and T0^2("°^e) from the intervals given in (g^l). The 
system of two equations is solved and the values of jr/j^p > and \rjii\'^ > thus obtained were used to 
obtained predictions for the half-life of the third isotope. The results for NMEs corresponding to the case 
"CD-Bonn, large, gA = 1.25" (see Table |T|) are given in Table HIl We note that the experimental lower 
bounds quoted in eq. (|4.3p have to be taken into account since they can further constrain the range of 
allowed values of and |7?_Rp. Indeed, an inspection of the values in Table HIl shows that not all the 
ranges predicted for the third half-life using the solutions obtained for |ryi{p and \ri^\'^ are compatible 
with the lower bounds on the half-live of the considered nuclear isotopes, given in (|4.3p . In this case, 



Table II: The predictions for the half-hfe of a third nucleus {A3, using as input in the system of equations for 
l?;^!^ and eq. (15. ip . the half-lives of two other nuclei {Ai, Zi) and {A2, Z-i). The three nuclei used are ^''Ge, 

^°°Mo and ^^"Te. The results shown are obtained for a fixed value of the half-life of (Ai, Z\) and assuming the 
half-life of (^2, Z2) to lie in a certain specific interval. The physical solutions for \r]v'^ and |77hP and then used 
to derive predictions for the half- life of the third nucleus (A3, Zz)- The latter are compared with the lower limits 
given in eq. (|4.3|) . The results quoted are obtained for NMEs given in the columns "CD-Bonn, large, qa = 1.25" 
in Table in One star beside the isotope pair whose half- lives are used as input for the system of equations (|5.ip . 
indicates predicted ranges of half- lives of the nucleus (^13,^3) that are not compatible with the lower bounds 
given in 
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some or all "solution" values of |77_rP and/or \rjy\'^ are ruled out. In Table HIl these cases are marked by a 
star. 

The results reported in Table |IT] are stable with respect to variations of the NMEs. If we use the 
NMEs corresponding to the case "CD-Bonn, large, qa — 1.0" (see Table U]), the limits of the 
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Figure 3: The values of the rescaled parameters \rj^'\^ (solid line) and |j7fl|^ (dashed lined), obtained as solutions 
of eq. (|5.ip for two values of T'i/2('^^G'e) and values of T'lj2C^^^o) lying in a specific interval. The physical 
(positive) solutions are delimited by the two vertical lines. The lower bound on Ty2C^^^Te) given in (|4.3|l does 
not lead to further constraints on and \riR\^ . The horizontal dashed line corresponds to the prospective 

upper limit from the upcoming H /3-decay experiment KATRIN [45l |. See text for further details. 



intervals quoted in Table HTl change by ±5%. If instead we use the NMEs corresponding to the Argonne 
potential, "large basis" and gA = 1-25 {qa = 1-0), the indicated limits change by ±10% (±14%). 

These results and considerations are illustrated in Figs. I2E1 The horizontal dashed line in these 
figures corresponds to the prospective limit planned to be obtained in the upcoming KATRIN experiment 
|45| . In figure [3] we show the solutions for |ryi{p and/or (conveniently rescaled), obtained for two 
values of T^J^{^^Ge) = 2.23 x lO^^ y and 10^^ y, assuming T^J^l}°°Mo) has a value in a certain interval. 
In the case of T^J^{^^Ge) = 2.23 x lO^^^ and T°J^(^^Ge) = lO^^ y, the derived physical values of \-qB? 
and |?7i,p lead to predictions for T^J^i^^'^Te) which are compatible with the existing lower limit (Fig. [31 
left and right panel). We get similar results using as input in the system of two equations for jTy^P 
and jr^yp the half- lives of different pairs of isotopes, and the lower limit of the half- life of the third as 




10- 



3.7x10^" 4.1x10^'' 4.5 xlO^"* 
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Figure 4: The same as in Fig. [S] but using as input hypothetical values of the half-lives of ^^Ge and ^'"'Te, 
T^j2C^Ge) and Tij2C^'^Te). The physical (positive) solutions are delimited by the two vertical lines. The lower 
bound on T{'y'2(^°°Afo) given in l|4.3p does not lead to further constraints on |77^,j;p. 
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ri/2("°Te) [y] r,/2(™Te) [y] 

Figure 5: The same as in Fig. [31 but using as input two values of the half-life of ^'"^Mo and values of the half-life 
of ^^"Te lying in a specific interval. The physical (positive) solutions are delimited by the two vertical lines. The 
lower bound on T^Y2C^Ge) given in (|4.3p does not lead to further constraints on [ry^P and |7)i?.p. 



an additional constraint. They are presented in Figs. In Fig. [7] we show the solutions for l?]^^ 

and \t]r\^ for T^J^{^^°Te) = 3 x 10^'' y and T°^^{^^°Mo) = (2.42 - 2.63) x lO^" y. In contrast to the 

cases illustrated in Figs. [3]- El most of the solution values of |?7^p and |r/;jp are excluded by taking into 
account the lower bound on Ty^C^Ge) given in eq. (|4.3p . 

We have studied also the dependence of the results discussed above on the value of gA and the NMEs 
used. This was done using the "large basis" NMEs, obtained with the CD-Bonn and Argonne potentials 
for the two values of gA — 1.25; 1.0. Some of the results of this study are presented graphically in Figs. 
[51 and ini The horizontal dashed line in these two figures corresponds again to the prospective limit from 
the upcoming KATRIN experiment jl^. We note that in the cases studied by us, changing the value of 
gA from 1.25 to 1.0 for a given potential (CD-Bonn or Argonne) 

does not lead to a significant change of the solutions for \r]^\'^ and the change is smaller than 

approximately 10%. The solutions exhibit a larger variation when for a given gA and basis, the NMEs 
calculated with the CD-Bonn potential are replaced by the NME's obtained with the Argonne potential 
(Fig. IHl upper right and lower left panels). In this case, as we have mentioned earlier, given Ti, the 
interval of allowed values of the half- life of the second nucleus T2 changes somewhat. In the specific 
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Figure 6: The same as in Fig. \5\ but using as input two values of the half-life of ^"^"Te and values of the half-life 
of ^""Mo lying in a specific interval. The physical (positive) solutions are delimited by the two thick vertical 
lines. The lower bound on r{%(^^G'e) does not lead to further constraints on It],,]^ and 
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Figure 7: The values of the rescaled parameters It^^P (solid line) and |j7fl|^ (dashed lined), obtained as solutions 
of eq. (|5.ip for the minimum value of T^j2C'^°Te) specified in eq. (14. 3|) . The physical (positive) solutions are 
delimited by the two vertical lines. The gray region is an excluded due to the lower bound on Tij2C^Ge) quoted 
in ((23]). 



cases shown, e.g., in Fig. [5] (upper right and lower left panels), Ti = T^^^C^Ge) = 2.23 x 10^^ y and the 
interval of interest of values of T2 = T^J^{^^'^Mo) shifts to larger values. Obviously, the solution values of 
the parameters jr^^p and |?7j^p, obtained with the two different sets of the NMEs, can differ drastically 
in the vicinity of the maximum and minimum values of T2, as is also seen in Figs. |5]and[ni If a given 
extreme value of T2, say max(T2), obtained with one set of NMEs, belongs to the interval of allowed 
values of T2 , found with a second set of NMEs, one of the fundamental parameters, calculated at max(r2) 
with the first set of NMEs can be zero, and can have a relatively large nonzero value when calculated 
with the second set of NMEs. Moreover, there are narrow intervals of values of T2 for which there exist 
physical solutions for jryi^p and |?7flp if one uses the NMEs obtained with the CD-Bonn potential and 
there are no physical solutions for the NMEs derived with the Argonne potential. If the measured value 
of T2 falls in such an interval, this can imply that either the two mechanisms considered are not at work 
in (/3/?)o,y-decay, or one of the two sets of NMEs does not describe correctly the nuclear transitions. As 
Figs. [5] and [5] indicate, the data from the KATRIN experiment can help limit further the solutions for 
|77i/p, obtained with NMEs calculated with the CD-Bonn potential and = 1.0 or with the Argonne 
potential. 

Let us note finally that Figs. |8]and|9]were obtained using hypothetical half-lives of '^^Ge and ^^°Mo. 



We get similar results if we use as input hypothetical half- lives of a different pair of nuclei, ^Ge and 
i30Te, i^o^e and ^°"Mo, etc. 
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Figure 8: Solutions for \ri^f (black lines) and \riRf (red lines), obtained by fixing Ti = T^J^C'^Ge) = 2.23 x 10^^ yr 
and Ta = T^}'2C°°Mo) and using the sets of NMEs calculated using the "large basis" and i) CD-Bonn potential, 
qa = 1.25 (solid lines) and qa = 1 (dashed lines) (upper left panel); ii) CD-Bonn (solid lines) and Argonne 
(dashed lines) potentials with qa = 1-25 (upper right panel); iii) CD-Bonn (solid lines) and Argonne (dashed 
lines) potentials with qa = 1-0 (lower left panel); iv) Argonne potential with qa = 1-25 (solid lines) and qa ~ 1 
(dashed lines) (lower right panel). The physical (positive) solutions for |?7,/p and Irmf shown with solid (dashed) 
lines are delimited by two vertical solid (dashed) lines. The horizontal dashed line corresponds to the prospective 
upper limit [isj |<m>| < 0.2 eV. 



VI. TWO "INTERFERING" MECHANISMS 



Neutrinoless double beta decay can be triggered by two competitive mechanisms whose interference 
contribution to the (/3/3)oiy-decay rates is non-negligible. In this Section we analyze the case of light 
Majorana neutrino exchange and gluino exchange. From equation (|4.2p it is possible to extract the 
values of |77i/p, |?7a'P and cos a setting up a system of three equation with these three unknowns using 
as input the "data" on the half-lives of three different nuclei. The solutions are given by: 

M^^^, \VX'\^ = ^, z = 2cosa|?7,||77A'| = (6.1) 
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Figure 9: The same as in Fig. HI but for Ti = T'^'J.^C'Ge) = 2.23 x lO^'^ yr and T2 = T°)'2C^°Te). 



where D, Di, D2 and Dt, are the following determinants 
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2, A' 
Of 

3. A' 



M' 

M' 



3,v 



{M'tf 



1/TiGi 
I/T2G2 
l/TaGs 



(6.3) 



We must require that 
ence term satisfies: 



\rji,\'^ and |?7a'P be non-negative and that the factor 2 cosa|77^||77A' | in the interfer- 



- < 2cosa|77,.||77A'| < 2\T]i,\\r]x'\. (6.4) 

If we fix (i.e. have data on) the half-lives of two of the nuclei and combine these with the condition in 
eq. (|6.4p . we can obtain the interval of values of the half- life of the third nucleus, which is compatible 
with the data on the half-lives of the two other nuclei and the mechanisms considered. The minimal 
(maximal) value of this interval of half-lives of the third nucleus is obtained for cos a = -1-1 (cos a = — 1). 
Examples of the intervals of half-life values of the third nucleus obtained using the half-life values of 



the other two nuclei ^ for the (/3/3)o;y-decay mechanisms discussed are Usted in Table IIIII The results 
reported in Table Hill are obtained with NMEs corresponding to the CD-Bonn potential, the "large basis" 
and gA = 1-25. 

Table III: Ranges of half-lives of T3 in the case of two interfering mechanisms: the light Majorana neutrino 
exchange and gluino exchange dominance. 



T'j";^ [y] (fixed) 


T^;;2 [y] (fixed) 


Allowed 


T(Ge) 


= 2.23 ■ 10^^ 


T(Afo) = 


= 5.8 ■ 1024 


5.99 


1024 < T{Te) < 7.35 ■ 


1024 


T(Ge) 


= 2.23 ■ 10^^ 


T(Te) = 


3 ■ 1024 


2.46 


1024 < T(Mo) < 2.82 


1024 


T(Ge) 


= 1026 


T(Mo) = 


= 5.8 ■ 1024 


6.30 


1024 < T{Te) < 6.94 • 


1024 


T(Ge) 


= 1026 


T(Te) = 


3 • 1024 


2.55 


1024 < T(Mo) < 2.72 


1024 


T(Ge) 


= 2.23 ■ 1025 


T(Te) = 


3 ■ 1025 


2.14 


1025 < T{Mo) < 3.31 


1025 


T(Ge) 


= 1026 


T(Te) = 


3 ■ 1025 


2.38 


1025 < T(Mo) < 2.92 


1025 



We show in few illustrative figures (Figs. fTUl - fO)) the results of the determination of \ri^\'^, |??a'P 
and cos a using different values of half-lives of the three nuclei ^^Ge, ^"^Mo and ^^°Te from the intervals 
given in eq. (|4.3p . The lower bounds of the half- lives quoted in eq. (|4.3p are taken into account. In these 
figures the physical allowed regions correspond to the areas shown in white, while the areas shown in 
gray are excluded. The allowed interval of values of the half-life of the 3rd nucleus, corresponding to the 
white areas, are listed in the 3rd column of Table Hill The results presented in Figs. [10] -[121 are derived 
using the NMEs, calculated with the CD-Bonn potential, the "large basis" and gA = 1.25. 

It is interesting to note that for the two fixed half- life values used to obtain Figs. [TUl [TTl and [T^ 
the interference between the contributions of the two mechanisms considered is destructive: one finds 
using these values that for most of the physical (positive) solutions for \'ri^\'^ and |??a'P, cos a is negative. 
Moreover, the rescaled parameters jT^i^p x 10^° and [t^a' P x 10^** in most of the solution regions have very 
close values. This is due to the fact that for most of the physically allowed values of \ri^\'^ and |?7a'P, 
each of the two terms including {rjul"^ or |?7a'P as a factor in the right hand side of eq. (|4.2p is much 
larger than the free term in the left hand side of eq. (j4.2p . As a consequence, in order to explain the 
"data" (i.e. the chosen values of the half-lives of the three isotopes) there should be a strong mutual 
compensation between the contributions due to the two mechanisms. This is possible only if |?7iyp(Af'j'j^)^ 

and |?yA' A')^ have close values and cos a = —1. 

In the case of destructive interference between the two contributions, |<to>| can have values which 
exceed the limit on the absolute scale of neutrino masses set by the ■^H /3-decay experiments |3|4^, eq. 
(|4.4p . This limit is indicated as a horizontal solid line in Figs. [TUl -[HI It leads to further constraints on 
the physical solution for |r/^P, and thus for |?/a'P- 

As we have already indicated, a more stringent limit on the absolute neutrino mass scale and therefore 
on I < TO > I is planned to be obtained in the KATRIN experiment fis']: it is given in eq. (|4.5p . The 
KATRIN prospective upper bound is shown as a horizontal dashed line in Figs. [TUl-lHl As the results 
presented in Figs. [TUl - [T^ indicate . if the limit of 0.2 eV will be reached in KATRIN experiment, if will 
lead to severe constraints on the solutions for \r]^\'^ obtained in the cases we have considered, strongly 
disfavoring (if not ruling out) essentially all of them. 

In Fig. [T21 we illustrate the possibility of constructive interference between the light neutrino and 
the gluino exchange contributions. The solutions shown in Fig. [121 are not constrained by the limits 
obtained in the '^H /3-decay experiments 0, HH; they also satisfy the prospective upper limit from 
KATRIN experiment. 

It is not difficult to derive from eqs. ()6.1[) - (|6.3p the general conditions under which |r/,/p and 
It^A'P are positive and the interference between the light neutrino and gluino exchange contributions is 
constructive (destructive), i.e. cos a (or z) is positive (negative). We will illustrate them later using 
again the NMEs calculated with the CD-Bonn potential, the "large basis" and gA = 1.25. 



5 Technically this is done in the following way. Fixing the half-lives of two isotopes, Ti and T2, and varying the half-life of 
the third isotope T3 in a certain interval, we obtained |»)j/|2, |r)_),/|2 and z = 2 cos a\ri,^\\rix' \ as a function of T3. Requiring 
that |»?i^|2 > 0, I^Ja'P > and that — 277,^1 |r;;^/ 1 < z < 2?7i^ 1 1)7;^/ 1 determines the interval of physically allowed values of T3 
(given Ti, T2 and the mechanisms of (/3/3)oi/-decay considered). This interval of physically allowed values of T3 is shown 
in Table ITnl 
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Figure 10: The values of the rescaled parameters |77^| (thick line) and \Tjx' \ (dashed lined), obtained as solutions 
of the system of equations ((42]) for fixed values of T^'j^C^Ge) and T^J^C^^Mo) and values of T°)'2C'^°Te) lying in 
a specific interval. The physical allowed regions correspond to the areas shown in white, while the areas shown 
in gray are excluded. The horizontal solid (dashed) line corresponds to the upper limit |<m>| < 2.3 eV 

(prospective upper limit [45] |<m>| < 0.2 eV). See text for details. 
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Figure 11: The same as in Fig. [TO]for a different set of values of the three half-lives used as input in the analysis. 
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Figure 12: The same as in Figs. [TU] and [TT] for a different set of values of the three half-lives used as input in the 
analysis. 



Consider first the conditions for constructive interference. We will introduce a somewhat simplified 
notations in this part of the article: Ti, Gi, Mi and for i = 1,2,3 will denote respectively the half-life, 
phase space factor, light neutrino and dominant gluino exchange NMEs for "^^Ge (i = 1), ^""Mo {i = 2) 
and i^OXe = 3). The first thing to notice is that it follows from Table|T]that the ratios of NMEs Af,;/A,; 




Figure 13: Left Panel: the values of 1 77^ ^ x 10^° (thick line), iTy^'pxlO^* (dashed line) and 2 = 2 cos a|77^| |?7a' I 10^^ 
(dot-dashed line) corresponding to the half-lives of ^^Ge, ^"''Mo and ^^"Te indicated on the figure. The interval 
of values of T^j2C^''Te) between the two vertical lines corresponds to physical (positive) solutions for 
and |77a'P a-s well as to a positive 2 (i.e. to a constructive interference between the contributions due to the 
two mechanisms). The horizontal solid (dashed) line corresponds to the upper limit (43 | (45l | |<m>| <2.3 eV 
(prospective upper limit [Zs^l |<7ti>| < 0.2 eV). See text for details. Right Panel: the phase a. 



satisfy the inequalities: 



Ml M2 M3 
Ai A2 A3 



(6.5) 



This inipHes that the determinant D, defined in eq. 



is negative: 



i. = A?A^A^(^-^)(^-^)(^-^), 
1 2 A2 Ai A3 Ai A3 A2 " 



(6.6) 



Consequently, in order to have \r]^\'^ >, |?7a'P > and constructive interference between the two contri- 
butions, i.e. z > 0, all three determinants Di, D2 and D^, defined in eqs. (|6.2p and (|6.3p . have to be 
negative: Da < 0, a = 1,2,3. Given the half- life Ti and the NMEs Mi and A^, these three conditions 
are satisfied if each of the two half-lives T2 and T3 lies in specific intervals ^: 



A) 



7^ Ti < < 



MiAi 
M2A2 
£2 



r;„ 1 



^MfA§H^gf(IifAi 



(M2A2A?-AfiAiAi) g|T2 



(M3 A3 A2 - Ml Ai A2 ) + (M2 A2 A2 - A/f3 A3 A| ) 



B) 



( MiAi Gl 
M2A2 G2 



Ti<T2<^^T,^ 



iM^Al-M^Al)^T2 



(MiAf-MfAl)^^ 



{MiAi 



{M2A2MI-M1A1MI) g| T2 



( A/f 3 A3 M 2 - Ml A 1 M| ) + ^gi. ( Af 2 A2 Af| - Af3 A3 Af| 



(6.7) 



For the NMEs calculated with the CD-Bonn potential, the "large basis" and gA — 1.25 and given Ti ^ 0, 



The quoted solutions are valid, as can be shown, provided M3/A3 < M2/A2 < 0.5(1 + v'SjAfs/Aa, which is fulfilled for 
the NMEs given in Table [l] 



the conditions for constructive interference, z > 0, 



z>0 : 



0.14 Ti < Ta < 0.16 Ti . 
0.16 Ti <T2 < O.lSTi, 



read: 
4.44 Ti Ta 



3.74 Ti - 0.93 
4.44 T1T2 



<T.< 



<T.< 



2.WT1T2 

1.78 Ti - 0.47 T2 
4.10 TiTa 

3.44 Ti - 0.81 T2 



(6.5 



3.74 Ti - 0.93 

These conditions imply that given Ti . a constructive interference is possible only if T2 lies in a relatively 
narrow interval and T3 has a value in extremely narrow intervals, the interval for T2 being determined 
by the value of Ti, while that for T3 - by Ti and the interval for T2. The fact that both the intervals 
for T2 and T3 are so narrow is a consequence of the values of the NMEs used, more precisely, of the 
fact that, for each of the two mechanisms discussed, the NMEs for the three nuclei considered differ 



relatively little: we have \Mi - Mj\ « Mi,Mj, \Ai - Aj\ « Ai,Aj, i ^ j = 1,2,3, and typically 
\Mi - Mj\/{0.5{Mi + Mj)) - 10~\ \Ai - Aj\/{0.5{Ai + A^)) - (lO'^ - IQ-^). We get similar results 
for the other sets of NMEs, quoted in Table |T] and calculated with the "large basis". In order to have a 
constructive interference in a much wider interval of values of T2, i.e., to have the minimal value of T2 
much smaller than the maximal value of T2 in case A) in eq. (j6.7p . for instance, the following inequality 
has to be satisfied: A1/A2 << M1/M2. An inspection of Table|T]shows that this inequality is not satisfied 
by any of the relevant sets of NMEs. Numerically, the intervals of values of T2 and T3 given in eq. (|6.8p . 
for which z > 0, are very similar to those quoted in eq. (j5.4p . 

For the value of T(^^Ge)= 2.23-10^^ y, for instance, the conditions for a constructive interference 
are given by: 



3.18 • 10^4 y < T2 < 3.55 • 10^* y 



24 



i.i9r2 



3.55 • lO^-^ y < T2 < 3.97 • 10^-* y 



24 



1.00- 1.12- 10-26 

1.186T2 



1T2 



1.00- 1.117- io-26y-ir2 



<r3 < 



1.19T2 



1.00- 1.18- 10-26y- 

i.i89r2 

1.00 - 1.059- 10 - 26y-lr2 



(6.9) 

Given the fact that 3.18 - 10^^ y < T2 < 3.97 - lO^"' y and that T2 enters in the denominators of the 



limiting values of T3 multiplied by 10 



the interval of values of of interest is extremely narrow. 



We have z > for, e.g., TC^Ge)= 2.23-10^^ y, T(i'^"Mo)= 3.7-10^^ y and T(""Te)= 4.58-10^* y, as is 
also illustrated in Fig. [T51 

There are cases in which one has 



\VX' 



can be derived from eqs. 



= or |77a'P = 0. The general conditions for having |77i^| 
([STT]) - dOJ and read: 

- Gi^ 7^, 







\VX' 



0,|r;,|VO 




(M2A2At-MiAiA^)g|T2 ^ 

(Af3A3A2-MiAiA|) + ^ga(A./2A2A§-Af3A3A|) ' 



(Af2A2Mi^-MiAiM|) T2 



( Af3 A 3 A/ 2 - A/i A 1 Af I ) + 4^ ( Af 2 A 2 A/ 1 - Af3 A 3 Af | ) 



(6.10) 



They correspond to some of the limiting values of T2 and T3 in eq. (|6.7p . We will illustrate them below 
numerically using the NMEs calculated with the CD-Bonn potential, the "large basis" and = 1.25. 
If, for instance, one fixes Ti = T(^^Ge)= 2.23-10^^, we have i) jr/i^p = (and zero interference term) for 
T2 = 3.18 - 10^4 y and T3 = 3.91 - 10^* y; ii) \r]y\^ = (and zero interference term) for T2 = 3.97 - 10^^ y 
and = 4.93 - 10^^ y, where T2 and T3 denote the half-lives of ^°°Mo and ^'^"Te, respectively. In general, 
given Ti we have \r]^\^ = 0, |?7a'P 7^ if 

2 10 T T 

T2 = 0.14T„ ^3= ^,3^^_;;,^^ ^0-18T,, (6.11) 
and |r7A'P = 0, \ri^\'^ ^ provided 

4.10 TiT-. 



The conditions for having zero inteference term, z — 0, but \rj^\^ 7^ or It^a'P 7^ 0, read: 



Ul rp ^ rp rp 



Mi 



(MiAi 



(A/|A?-M2Ai) + if^(M|Ai-M|Ai) 



(6.13) 



Given Ti, the general conditions for destructive interference, i.e. for z < 0, can be derived in a similar 
way. They read: 



A) 



B) 



< T2 < 



o<n< 



(M2A2A^-J\/iAiA^) ■ 



■T2 



(A/aAsAf-MiAiA: 



M3A3AI) 



A? 



< T, < 



0<T, < 



(M^^A' 



(M|A2, 



-M^Ai)+"^ 



(M|Ai 



(6.14) 



.15) 



C) 



Gl rp rp M1A1M3-A3MI Gl rp 



G 

0<T3 < 



(Ma A2 Mj^ - Ml Ai M| ) 



G2 



T2 



(M3A3M=-MiAiM|) + ^ 



(M2A2Af|-M3A3Af|) ' 



(6.16) 



- A/2A2Af3-A3A'/| G2 1 ■ 



T3 >0 



(6.17) 



Obviously, one has to take into account the existing experimental lower limits on T2 and in eqs. (|6.14p 
- (|6.17p . We will give next the "numerical" equivalent of the conditions (|6.14p - (|6.17p . obtained with 
NMEs calculated with the CD-Bonn potential, the "large basis" and gA = 1.25: 



z < : <^ 



T2<0.UTi, T:,< 

0.14 Ti <T2< 0.18 Ti 

0.18 Ti < T2 < 4.23 Ti 
T2>4.23ri r3>o. 



2.IOT1T2 



1.78Ti - O.47T2 ' 
4.44 Ti T2 

3.74Ti -O.93T2 
4.10 riT2 

3.44 Ti -0.81 T2 



(6.18) 



The intervals of values of T2 and T3 in eqs. (|6.14p - (|6.18p are very different from those corresponding to 
the case of two "non-interfering" (/3/3)o,y-decay mechanisms given in eq. (j5.4p . the only exception being 
the second set of intervals in eq. ()6.18|) , which partially overlap with those in eq. ()5.4p . This difference 
can allow to discriminate experimentally between the two possibilities of (/3/3)ojy-decay being triggered 
by two "non-interfering" mechanisms or by two "destructively interfering" mechanisms. We have check 
how the intervals of values of the half- life T3 given in Table IIIII corresponding to NMEs derived with the 
CD-Bonn potential, gA — 1-25 and the "large basis", change when one uses the NMEs obtained with the 
same potential and basis, but using gA — 1-0, as well as the NMEs found with the Argonne potential for 
gA — 1.25; 1.0 and the "large basis". The results are shown in Tables HVl - IVII We see that for certain 
values of the hypothetical half-lives of the two nuclei, the interval of allowed values of the half- life of 
the third nucleus becomes noticeably larger when calculated with NMEs, corresponding to gA = 1.0 
or to the Argonne potential. This is due to a relatively deep compensation between the three terms 
in the (/3/3)o;/-decay rate of the third nucleus in the case of a negative interference term (destructive 
interference) . 



Similar analysis can be performed for any other pair of "interfering" mechanisms assumed to be 
operative in (/?/?) ojy-decay. We note also that the extension of the analysis to more than two mechanisms 
generating the (/3/3)oi/-decay is rather straightforward. 



Table IV; CD-Bonn potential and qa = 1 



W (fixed) 


T^^a (fixed) 


Allowed 


T(Ge) = 2.23 • 10^5 
T(Ge) = 2.23 • 1025 
T(Ge) = 2.23 • 1025 
T(Ge) = 10^6 
T(Ge) = 10^6 
T(Ge) = 1026 


T(Mo) = 5.8 • 10^4 
T(Te) = 3 • 1024 
T(Te) = 3 • 1025 
T(Mo) = 5.8 • 1024 
T(Te) = 3 • 1024 
T(Te) = 3 • 1025 


3 • 1024 < T(Te) < 8.62 • 1024 
2.55 • 1024 < T{Mo) < 6.18 • 1024 
1.33 • 1025 < T(Mo) < 3.88 • 1026 
3.62 • 1024 < T{Te) < 6.04 • 1024 
3.11 • 1024 < T{Mo) < 4.70 • 1024 
2.15 • 1026 < T(Mo) < 8.29 • 1025 


Table V: Argonne potential and qa = 1-25 


T^;;2 (fixed) 


tl'j^ [y] (fixed) 


Allowed 


T(Ge) = 2.23 • 10^5 
T(Ge) = 2.23 • 1025 
T(Ge) = 2.23 • 1025 
T(Ge) = 1026 
T(Ge) = 1026 
T(Ge) = 1026 


T(Mo) = 5.8 • 10^* 
T(Te) = 3 • 1024 
T(Te) = 3 • 1025 
T(Mo) = 5.8 • 1024 
T(re) = 3 • 1024 
T(Te) = 3 • 1025 


3 • 10^* < T(Te) < 9.22 • 10^* 
2.55 • 1024 < T(Mo) < 7.92 • 1024 
1.19 • 1025 < T(Mo) < 2.55 • 102'^ 
3.15 • 1024 < T{Te) < 5.85 • 1024 
3.25 • 1024 < T(Mo) < 5.49 • 1024 
2.08 • 1025 < T{Mo) < 1.20 • 1026 



Table VI: Argonne Potential and pA = 1 



T'ji^j [y] (fixed) 


T^;;2 (fixed) 


Allowed 


T(Ge) = 2.23 • 10^5 
T(Ge) = 2.23 • 1025 
T(Ge) = 2.23 ■ 1025 
T(Ge) = 102e 
T(Ge) = 1026 
T(Ge) = 1026 


T(Mo) = 5.8 • 10^" 
T(Te) = 3 • 1024 
T(Te) = 3 • 1025 
T(Mo) = 5.8 • 1024 
T(Te) = 3 • 1024 
T(Te) = 3 • 1025 


3 • 10^* < T(Te) < 1.11 • 10^5 
2.63 • 1024 < T{Mo) < 2.04 • 1025 
9.19 • 1024 < T{Mo) < 2.36 • 1026 

3 • 1024 < T(Te) < 5.07 ■ 1024 
3.82 • 1024 < T(Mo) < 9.44 • 1024 
1.96 • 1025 < T{Mo) < 6.54 • 1026 



Finally, we would like to point out to one additional consequence of the "positivity" conditions and 
the condition the interference term should satisfy when two "interfering" mechanisms are responsible for 
the (/3/3)ojy-decay. Let us denote the two fundamental parameters characterising the two mechanisms 
by r]p and 77^. Then, given the half-life of one isotope, say of ^^Ge (Ti), and an experimental lower 
bound on the half-life of a second isotope, e.g., of ^^'^Te (Ts), the conditions |?7^P > 0, jr^^P > and 
~I'7/3||^k| — l^/sll^yftl cosa^K < |77^||?7k|, imply a constraint on the half-life of any third isotope, say of 
^°°Mo (T2). This latter constraint depends noticeably on the type of the two "interfering" mechanisms 
generating the (/3/3)oi/-decay and can be used, in principle, to discriminate between the different possible 
pairs of "interfering" mechanisms. Below we illustrate this result by deriving the constraint one obtains 
on the half-life of ^°°Mo, T2, assuming that the half-Hfe of ^^Ge is Ti = 2.23 x lO^^ y and taking into 
account the current experimental lower bound on the half-life of ^'^"Te, T3 > 3.0 x 10^'' y. Using these 
"data" as input, the NMEs calculated with the CD-Bonn and Argonne potentials, the "large basis" and 
gA — 1.25, we get the following constraint on T2 for the different pairs of "interfering" mechanisms 
discussed by us (the numbers in brackets are obtained with the NMEs corresponding to the Argonne 
potential, unless otherwise indicated). 

Light Neutrino and gluino exchange mechanisms: 

T2 = T^)'2C°°Mo) > 2.46 (2.47) x 10^'' y. (6.19) 

Increasing the value of T^^^C^Ge) leads to the increasing of the value of the lower limit. 

Light Neutrino and LH Heavy neutrino exchange mechanisms: 

T°;2(^°"Mo) > 2.78 (2.68) x lO^" y. (6.20) 

The value of the lower limit increases with the increasing of the value of the half-life of ^^Ge. 

LH Heavy neutrino and gluino exchange mechanisms: 

1.36 X 10^^ y < T°/2(^™Mo) < 3.42 x 10^'' y . (6.21) 

Increasing the value of T^^^C^Ge) leads to a shift of the interval to larger values, and for a sufficiently 
large T°^^C^Ge) > 10^6 y - even just to a lower bound on T°J^(^°^^Mo). For Ti = lO^^ y, for instance, 
we find 4.19 x lO^^ y < T^^^i^^^Mo) < 3.39 x lO^^ y. 

Using the NMEs derived with the Argonne potential we find a very different result - only a lower 
bound: T^J^C^^^Mo) > 5.97 x 10^^ y. The difference between the results obtained with the two sets of 
NMEs can be traced to fact that the determinant D in eqs. (|6.ip and (|6.2p . calculated with the second 
set of NMEs, has opposite sign to that, calculated with the first set of NMEs. As a consequence, the 
dependence of the physical solutions for jjy^P and on Ti, T2 and T3 in the two cases of NMEs is 

very different. 

Squarks-neutrino and gluino exchange mechanisms: 

Ti/2C°°Mo) > 7.92 (22.1) x 10^^ y. (6.22) 

For larger values of T^^^(^^Ge), this lower bound assumes larger values. 

We see that the two sets of NMEs lead to quite different results in the cases of the LH heavy 
neutrino and gluino exchange and squarks-neutrino and gluino exchange mechanisms. Nevertheless, the 
constraints thus obtained can be used, e.g., to exclude some of the possible cases of two "interfering" 
mechanisms inducing the (/3/3)oi/-decay. Indeed, if, for instance, it is confirmed that Tjy'jl^^Ge) = 
2.23 X 10^^ y, and in addition it is established, taking all relevant uncertainties into account, that 
T^j'^(^OQj^,j^^ ^q24 ^T^g,^ combined with the experimental lower limit on T^J^{^'^^Te) would rule out 

i) the light neutrino and gluino exchanges, and ii) the light neutrino and LH heavy neutrino exchanges, 
as possible mechanisms generating the (/3/3)oi/-decay. 



VII. SUMMARY AND CONCLUSIONS 



In the present article we have considered the possibility of several different mechanisms contributing 
to the (/3/3)oi/-decay amplitude in the general case of CP nonconservation. The mechanisms discussed are 
light Majorana neutrino exchange, exchange of heavy Majorana neutrinos coupled to (V-A) currents, ex- 
change of heavy Majorana neutrinos coupled to (V+A) currents, lepton charge non-conserving couplings 
in SUSY theories with _R-parity breaking. Of the latter we have concentrated on the so-called "domi- 
nant gluino exchange" mechanism. Each of these mechanisms is characterised by a specific fundamental 
lepton number violating (LNV) parameter. The latter are defined in Section 2. We have investigated in 
detail the cases of two "non-interfering" and two "interfering" mechanisms, generating the (/3/3)o;/-decay. 
In the analysis we have performed, we have used hypothetical (/3/3)oiy-decay half-lives of the following 
three isotopes: "^^Ge, '^'^"Mo and ^'^°Te. They are denoted as Ti, T2 and T3, respectively. Four sets of 
nuclear matrix elements (NMEs) of the decays of these three nuclei were utilised: they were obtained 
with two different nucleon-nucleon potentials (CD-Bonn and Argonne) and two different values of the 
axial coupling constant gA = 1.25; 1.0 (see Table |T1). 

If the (/3/3)oy-decay is induced by two "non-interfering" mechanisms, which for concreteness we have 
considered to be the light LH Majorana neutrino exchange and the heavy RH Majorana neutrino exchange 
with (V + A) currents, one can determine the squares of the absolute values of the two LNV parameters, 
characterising these mechanisms, Ir^j^p and \riu\'^, from data on the half-lives of two nuclear isotopes. 
We have done that using as input all three possible pairs of half-lives of ^^Ge, -'^""Mo and '^'^"Te, chosen 
from the intervals given in eq. (|4.3p and satisfying the existing experimental lower limits, as well as the 
half-life of the (/3/3)o^-decay of '^^Ge, claimed to be observed in [11]: T^J^C^Ge) = 2.231°;;^^ x lO^^ y. 

We find that if the half-life of one of the three nuclei is measured, the requirement that |?7iyp > and 
\i1r\^ ^ ("positivity condition") constrains the other two half-lives (and the (/3/3)oi/-decay half-life of any 
other (/3/3)oiy-decaying isotope for that matter) to lie in specific intervals, determined by the measured 
half- life and the relevant NMEs (see eqs. (|5.3p - (|5.6p ). This feature is common to all cases of two "non- 
interfering" mechanisms generating the (/3/3)oi/-decay. The indicated specific half-life intervals for the 
various isotopes, are stable with respect to the change of the NMEs (within the sets of NMEs considered 
by us) used to derive them. The intervals depend, in general, on the type of the two "non-interfering" 
mechanisms. However, these differences in the cases of the (/3/3)o,y-decay triggered by the exchange of 
heavy Majorana neutrinos coupled to (V+A) currents and i) the light Majorana neutrino exchange, 
or ii) the gluino exchange mechanism, or i ii) the squark-neutrino exchange mechanism, are extremely 
small. One of the consequences of this feature of the different pairs of "non-interfering" mechanisms 
considered by us is that if it will be possible to rule out one of them as the cause of (/3/3)oiy-decay, most 
likely one will be able to rule out all three of them. Using the indicated difference to get information 
about the specific pair of "non-interfering" mechanisms possibly operative in (/3/3)oi/-decay requires, in 
the cases considered by us, an extremely high precision in the measurement of the (/3/3)oiy-decay half-lives 
of the isotopes considered. The levels of precision required seem impossible to achieve in the foreseeable 
future. If it is experimentally established that any of the indicated intervals of half-lives lies outside the 
interval of physical solutions of \ri^\'^ and \t]r\'^, obtained taking into account all relevant uncertainties, 
one would be led to conclude that the (/3/3)oiy-decay is not generated by the two mechanisms considered. 
The constraints under discussion will not be valid, in general, if the (/3/3)oi/-decay is triggered by two 
"interfering" mechanisms with a non-negligible (destructive) interference term, or by more than two 
mechanisms none of which plays a subdominat role in (/3/3)o!y-decay. 

We have studied also the dependence of the physical solutions for jTy^p and |r7_Rp obtained on the 
NMEs used. Some of the results of this study are presented graphically in Figs. [5] and [51 We found that 
the solutions can exhibit a significant variation with the NMEs used. Given the half-life Ti, the interval 
of allowed values of the half-life of the second nucleus T2 , determined from the "positivity conditions", 
^ 0, |??_Rp > 0, changes somewhat with the change of the NMEs. The solution values of the 
parameters |r?yp and |?7i{p, obtained with the two different sets of the NMEs, can differ drastically in 
the vicinity of the maximum and minimum values of T2 (Figs. [5] and [5]). If a given extreme value of T2, 
say max(T'2), obtained with one set of NMEs, belongs to the interval of allowed values of T2, found with 
a second set of NMEs, one of the fundamental parameters, calculated at max(T2) with the first set of 
NMEs can be zero, and can have a relatively large nonzero value when calculated with the second set of 
NMEs. Moreover, there are narrow intervals of values of T2 for which there exist physical solutions for 
jryi/p and |?7ji;p if one uses the NMEs obtained with the CD-Bonn potential and there are no physical 



solutions for the NMEs derived with the Argonne potential. If the measured value of T2 falls in such an 
interval, this can imply that either the two mechanisms considered are not at work in (/3/3)o,y-decay, or 
one of the two sets of NMEs does not describe correctly the nuclear transitions. 

Neutrinoless double beta decay can be generated by two competitive mechanisms whose interference 
contribution to the (/3/3)oiy-decay rates is non-negligible. In the case when two "interfering" mechanisms 
are responsible for the (/3/3)oi/-decay, the squares of the absolute values of the two relevant parameters 
and the interference term parameter, which involves the cosine of an unknown relative phase of the 
two fundamental parameters, can be uniquely determined, in principle, from data on the half-lives of 
three nuclei. We have analyzed in detail the case of light Majorana neutrino exchange and gluino 
exchange. In this case the parameters which are determined from data on the half-lives are \rj^\'^ , ^A'P 
and z = 2 cos a\ri,j\\rixi\. The physical solutions for these parameters have to satisfy the conditions 
\ilv\'^ — 0; ^A'P > and — 2|77iy||?7A' | < z < 2\r],j\\r]\'\. The latter condition implies that given the 
half- lives of two isotopes, Ti and T2, the half-life of any third isotope T3 is constrained to lie is a specific 
interval, if the mechanisms considered are indeed generating the (/3/3)o,y-decay. If further the half-life of 
one isotope Ti is known, for the interference to be constructive (destructive) the half-lives of any other 
pair of isotopes T2 and T3, should belong to specific intervals. These intervals depend on whether the 
interference between the two contributions in the (/3/3)o,y-decay rate is constructive or destructive. We 
have derived in analytic form the general conditions for i) constructive interference (z > 0), ii) destructive 
interference (z < 0), iii) |r?^|2 = 0, It^a'P ^ 0, iv) jry^p ^ 0, |?7a'P = and v) z = 0, \r]^\^ 0, \r]y\^ ^ 0. 

We have found that, given Ti, a constructive interference is possible only if T2 lies in a relatively 
narrow interval and T3 has a value in extremely narrow intervals, the interval for T2 being determined 
by the value of Ti , while that for T3 - by Ti and the interval for T2 . The fact that both the intervals for 
T2 and T3 are so narrow is a consequence of the fact that, for each of the two mechanisms discussed, the 
NMEs for the three nuclei considered differ relatively little: the relative difference between the nuclear 
matrix elements of any two nuclei does not exceed 10%. 

The intervals of values of T2 and T3 corresponding to destructive interference ( eqs. (|6.14p - (|6.18p ) 
are very different from those corresponding to the cases of constructive interference and of the two 
"non-interfering" (/3/3)o,y-decay mechanisms we have considered fea. (j5.4|) ). Within the set of {f3/3)o,y- 
decay mechanisms studied by us, this difference can allow to discriminate experimentally between the 
possibilities of the (/3/3)o;y-decay being triggered by two " destructively interfering" mechanisms or by 
two "constructively interfering" or by two "non- interfering" mechanisms. 

We have shown also that further significant constraints on the physical solutions for the fundamental 
parameter \ri^\'^ in the case of the light Majorana neutrino exchange mechanism and the gluino exchange 
(or any other "interfering") mechanism can be obtained by using the current and the prospective upper 
bounds on the absolute scale of neutrino masses from the past (fl IZsj and the upcoming KATRIN 
(4^ /3-decay experiments of 2.3 eV and 0.2 eV, respectively. Our results show that the KATRIN 
prospective upper bound of 0.2 eV, if reached, could imply particularly stringent constraints in the cases 
of "destructively interfering" mechanisms one of which is the light neutrino exchange, to the point of 
strongly disfavoring (or even excluding) some of them. The KATRIN prospective upper bound could be 
used to constrain also the fundamental parameters of two "non-interfering" mechanisms, one of which 
is the light Majorana neutrino exchange. This bound could eliminate, in particular, some parts of the 
half- life solution intervals where there is a significant dependence of the values of |77^p obtained on the 
set of NMEs used. 

The measurements of the half-lives with rather high precision and the knowledge of the relevant 
nuclear matrix elements with relatively small uncertainties is crucial for establishing whether more than 
one mechanisms are operative in (/3;S)oi/-decay. The method considered by us can be generalised to the 
case of more than two (/3/3)o^-decay mechanisms. It allows to treat the cases of CP conserving and CP 
nonconserving couplings generating the (/?/3)oj,-decay in a unique way. 
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Note Added. The possibility of several mechanisms being operative in (/?/3)oiy-decay is also dis- 
cussed in another recent preprint |48i] , where the sets of nuclear matrix elements given in Table U are also 
used. However, the aspects of the problem of multiple mechanisms generating the (/3/3)o;y-decay investi- 



gated in the present article and in the preprint |48l] are very different and, apart from the description of 
the nuclear matrix elements, the two studies practically do not overlap. 
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